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ABSTRACT 


A study of the physics of the piano reveals that while the upper 
partials of the steel strings are the eigen-frequencies of the complex 
tone, they are not integer multinles of the resnective fundamentals. 
To properly measure and tune these eigen-partials, a digital filter 
capable of sweeping a major portion of the audio-frequency spectrum 
had to be implemented. Such a filter, a tuneable fixed-coefficient 
digital filter, is discussed as well as a simple pole-zero desiqn 
procedure for determining the required coefficients. Each module, 
including the Frequency Deviation Detector and Counter, the Time-Base 
Generator, the Digital Filter, the Reference Frequency Generator and 
the Display and Control Module, of the proposed tuner is illustrated 


and discussed, 
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I. INTRODUCTION 


The purpose and end result of this paper is to demonstrate how 
a digital piano tuner may be constructed. After only a small amount 
of investigation, however, it became readily apparent that a certain 
understanding of not only the basic physics of the piano, but also 
of music theory, was essential before a true definition of the 
problem could be established. 

It was relatively simple to find many books on music and sound 
that gave the precise frequencies for each note of the "Well-Temper- 
ed Piano"! (See Fig. 1). Using these frequencies it proved no 
great task to design and build a digital circuit that would in fact 
measure, with a very high degree of accuracy, these various notes .° 
It was soon found, however, that if this was done, such tuning would 
produce complete aural chaos to both the trained and untrained ear 
alike. As asserted above, to fully appreciate this "discovery" 


a small amount of musical and sound theory must be presented. 


A. SOUND AND PHYSICS OF THE PIANO 
When any string is struck or plucked and then allowed to vibrate, 
it not only vibrates at its fundamental frequency but also at its 


various upper harmonics or partials [Refs. MF-3, MJ-10 and MO-14]. 


Veach established the scale of equal temperament as an accepted 
scale in the musical world. His pieces for piano were all written in 
this scale. For a discussion of the history and explanation of the 
scale itself, see Refs. MH-5, M)-14, and MW-24. 


ese section head "A Frequency Deviation Detector and Counter Module." 


The fundamental frequency is determined by the length of the wire or 
string being considered and at what place along its length it is 
struck [Ref. MJ-10]. The upper partials occur because of a certain 
amount of stiffness in almost all strings, especially the steel 
strings used in the piano. 

Fig. 2 illustrates the vibration of an ideal string, that is, 
one without any stiffness. This string can be made to vibrate at 
many different frequencies. The fundamental frequency (a) produces 
a pure tone rarely heard in music. The higher-pitched partial tones, 
or overtones, are produced by harmonic vibrations (b) and (c), whose 
frequencies are integer multiples of the fundamental frequency. 

Fig. 3 illustrates the simultaneous vibration of a string at 
two or more different frequencies and illustrates the mode that is 
normal for all stringed instruments. 

The complex sound produced by this combination of separate tones 
has a timbre, or characteristic quality, that is determined largely 
by the number of partial tones and their relative loudness [Refs. MB-1, 
MC-2, MJ-10, MM-12, MP-19 and MO-14]. It is timbre that enables one 
to distinguish between two musical tones that have the same pitch 
and loudness, but are produced by two different musical instruments. 
With this in mind, it seems reasonable to conclude that the upper 
partials and not the fundamental frequency are the determining factors 
in the sound that is perceived [Ref. MC-2]. 

In order to further substantiate this conclusion, a test done by 
Bell Telephone at the Smithsonian Institute several years ago was 


investigated. The Bell engineers made two recordings of musical 


instruments, opera singer, noises and everyday sounds. In one, the 
fundamental frequency was extracted, while the other was left untouched. 
When played before both laymen and professionals in different fields, 
the groups were unable to determine which recording had the missing 
fundamental! In fact, even when all the frequencies of a musical 
composition below 300 hertz were removed, the quality of the music 
remained the same to a surprising degree [Ref. MC-4]. The "case of 
the missing fundamental" has been attributed to the fact that the 
ear uses periodicity and not frequency as a basis for pitch percep- 
tion [Ref. MP-18]. 

Another interesting fact of similar nature is brought out in 
Refs. MC-2, and MJ-10. They show, for the modern piano, that although 
the fundamental is the "loudest" when the key is first struck it 
rapidly dies out and upper partials, having a much longer decay time, 
take over the pitch of the note. (In very old pianos, the fundamen- 
is actually absent!) All of this further substantiates the conclu- 
sion above that the upper partials are the most important fre- 
quencies in a complex musical sound. 

"So what?," may very well be the question in the mind of the 
reader at this time. With a seemingly logical approach, it could 
very well be asserted that when the fundamental is measured, the 
upper partials are also measured and once again only a relatively 
simple frequency measuring scheme is required. This is, however, 
not the case. 

As stated before, piano strings are made of steel because of 


the tremendous stresses they are required to withstand in order to 


be able to generate an appreciable amount of sound. This use of 

a "wire" or steel string, unlike the gut string of the violin, 
introduces a certain amount of stiffness or inelasticity. This 
inelasticity has the very undesirable effect on the upper partials 
of causing them to vibrate at frequencies that are not integer 
multiples of the fundamental frequency. No longer can the formula: 


J = nf (1) 


be applied in determining the ro partial. The formula is shown 


in Ref. MF-3 to be: 


f = nf,((1+Bn*) s (7 +B) 2 = (2) 
where B is dependent on the dimensions of the wire! B is in fact 

a function of the length squared, diameter, Young's modulus of 
elasticity, the area of the cross section, the radius of gyration 
squared and the tension of the wire when jn the neutral positton. 

As can be easily imagined this causes the upper partials to differ 
from piano to piano and makes completely useless any attempt to tune 
by zero beating the fundamentals. > It should be noted here that the 
inharmonicity caused by this string stiffness always makes the 

upper partials sharp with respect to the desired, pure, integer 


harmonic. It is less pronounced in the lower octaves and greater 


jn the higher registers. 


Sin Ref. MM-13, Professor Franklin Miller suggests a means 
whereby the desirable result of reducing this inharmonicity of 
partials to a negligible quantity by applying a small amount of 
mechanical loading near one end of the piano string. 


|Z 


With the above comments made, it appears appropriate to make a 
small digression into the field of music theory and see what effect 


this inharmonicity has on a musical score. 


B. MUSIC THEORY 
The average human ear can distinguish about 1,400 discrete 
frequencies. However, in the equally tempered scale covering the 
hearing range from 16 to 16 x 10° hertz, there are only 120 dis- 
crete tones. 
Comparatively few people are capable of recognizing the true 
pitch of a musical tone. A great many individuals are able, though, 
to distinguish the ratio of two frequencies. Furthermore, most 
people recognize that when two notes are sounded together or immedi- 
ately following one another they either produce a pleasing effect or 
give a decidedly unpleasant reaction. In music theory these reactions 
are termed consonance and dissonance, respectively. The frequency 


ratios that produce these conditions are termed musical intervals or 


Simply intervals, 


+n though really beyond the scope of this paper, the equal tem- 
pered scale brought to the fore-front by Bach in his pieces for the 
“Well Tempered Klavier" is of great importance in music theory. The 
scale allows the musician or composer to modulate or change keys of a 
composition without changing the entire frequency spectrum of the 
piano. This latter change in the frequency spectrum for every change 
in key was why the more pleasing-to-the-ear scale of just intonation 
was finally dropped. It is interesting to note that all music before 
Bach's time was written in the scale of just intonation and modern 
day audiences are therefore not hearing these pieces as composed 
because of the change in scale. For a much deeper and interesting 
discussion of this subject see Refs. MH-6, MJ-9, MO-14, MP-15, 
Mw-21 and MW-24. 


(ls: 


The question now arises: why when two notes e.qg., C4 (261.63 Hz) 
and D4 (293.66 Hz) are sounded together a very rough sensation 1s 
produced, while when C4 and E4 (329.63 Hz) are played, a very smooth 
and pleasing musical effect occurs. The answer was found first by 
the great investigator Helmholtz: 

When two musical tones are sounded at the 

same time, their united sound is generally 

disturbed by the beats of the upper par- 

tials (harmonics), so that a greater or 

less part of the whole mass of sound is 

broken up into pulses of tone, and the joint 

effect is rough. This relation is called 

dissonance. But there are certain deter- 

minate ratios between frequency numbers, 

for which this rule suffers an exception, 

and either no beats at all are formed, or 

at least only such as have so little inten- 

sity that they produce no unpleasant dis- 

turbance of the united sound. These 

exceptional cases are called consonances. 
Helmholtz also found that this "roughness" was maximized at 33 beats/ 
sec and went as high as 132 beats/sec before its effect on musical 
sound became unappreciable. 

Now returning to the original question, concerning the differ- 
ence between C4-D4 (a musical second) and C4-E4 (a musical third), 
let the whole notes in Fig. 4 represent the notes played and the 
stars their respective upper partials. 

It becomes readily apparent from the illustration that the beating 
of upper partials of the second produce a beat frequency well within 


Helmholtz's "critical area," whereas the musical third does not. 


“Helmholtz .H.L.F. On the Sensations of Tone, 2d Eng. ed., 
trans. Ellis, A., p. 179-197, Dover Publications, 1954. 
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Fig. 4. Beats produced by a dissonance interval (a) 
and a consonance interval (b). 








The intervals that produce pleasing sensations have long been 
known to musicians and composers. It is the mixina of these various 
degrees of smooth harmonious sounds by themselves and with some of the 
dissonance intervals that creates all of the music heard today. Table 1. 
gives the various consonant intervals in order of decreasing agreeable 
audi tory a Notice the appearance of the musical third, while 


the second is conspicuously absent. 


NAME RATIO NAME RATIO 


Unison Lael Minor Third 62905 

Octave Co Major Sixth 58 

Fifth Se 2 Minor Sixth Giro 

Fourth C3 Major Seventh lSe6 
ogee! 


Major Third Minor Seventh 9:35 
Table 1. Pairs of tones producing agreeable auditory effects. 


Before attempting to tie all of this together, a recent investiga- 
tion by Plomb [Ref. MP-18] will lay the foundation for the binding. 
Plomb conducted this investigation to determine what factor or factors 
in a complex sound caused a listener to assign it a definite pitch. 

The results are therefore not only applicable to the complex sounds 
generated by the piano but can be applied to any complex tone. Plomb 
found that for fundamental frequencies up to about 350 hertz, the nitch 
was determined by the fourth and higher partials; for frequencies up to 
about 700 hertz, by the third and higher partials; and for frequencies 
up to about 1400 hertz, the complex tone was determined by the second 
and higher partials. In all of the above cases, the fundamental was 
shown to have no bearing on the determination of the pitch once the 
tone had been generated. For frequencies above 1400 hertz, however, 
the fundamental appeared to be the determining factor. This was 
attributed to the fact that the ear starts having trouble detecting 

the periodicity of the tone at these high frequences. 

It now takes little imagination to guess what would happen if the 
previously discussed inharmonicity of the piano wires were taken into 
account in the example of the second and third (see Fig. 4). The upper 
partials of D4 would need to be sharpened by only a very small? amount 
before this interval would cross Helmholtz's critical area and become 
quite rough and unpleasant to the listener. 

Thus the question arises: Just how much does the stiffness of the 


string cause the upper harmonics to vary? This question is easily 
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answered graphically in Fig. 5 [Ref. MB-1]. As is readily apparent, 


the higher partials can be as much as two full semi-tones sharp.° 


C. STATEMENT OF THE PROBLEM 

In the previous discussion it has been shown that while all musical 
strings vibrate at their fundamental frequency, they also vibrate at 
their various upper partials or harmonics. It was also pointed out 
that the upper partials and not the fundamental determines the pitch 
of the complex sound. This dependency on harmonics would present no 
real problem in the digital tuning of the piano if the strings of the 
piano vibrated at partials that were integer multiples of the fundamen- 
tal. This, however, was shown not to be the case. Therefore in order 
to properly tune the steel-stringed piano, one of the unper partials, 
as determined by Plomb, must be filtered from the complex sound, meas- 
ured and precisely tuned as if it exhibited the characteristics of a 
pure harmonic. 

1. A Commercial Analog Device 

Conn Instruments Inc. now has the only instrument on the market 

that attemots to solve this problem [Ref. MK-l11]. It is an analog 
calculator and filter that will measure the upner partials of the 
various keys of the piano. However, Conn admits that its tuner, "The 


Strobotuner," requires the aid of a good pair of musical ears.’ 


On semi-tone is the difference in frequency between two adjacent keys 
in the equal-tempered scale. This amounts to about a six per cent differ- 
ence in frequency between notes. Musicians call 1/100 of a semi-tone a 
Cenc? 


it is not desired to degrade from the quality of this instrument in 
any way. It is in fact a very fine machine that goes a long way in at- 
tempting to solve a very complex problem. Further information may be ob- 
tained by writing Conn Instruments Inc. Elkhart, Indiana 46514. 
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This author believes this is due to the fact that the bearing is laid by 
measuring the fundamentals of the octave chosen and then using nothing 
higher than the first partial in all subsequent tuning.© As Plomb's 
investigation has already shown, this is by no means sufficient for 
the majority of the notes on the piano. 

This writer decided that a better tuner could be built if one 
would filter, measure and tune the partials that Plomb's investigation 
proved to be of prime importance. 

2; The Author's Proposed Implementation 

The original goal was the desian of a digital piano tuner. 
This goal resulted in the decision that some sort of digital filter 
had to be implemented. Although there has been a great deal of work 
done in the field of digital filtering [Refs. EB-4, EC-7 thru EC-10, 
EG-15, EG-17 thru EK-19, EK-21, EM-28, EN-30, ER-33, ER-37, ER-38, 
ET-42 thru EV-46, EW-48 and EW-49], the area of tuneable fixed-coef- 
ficient digital filtering has been seemingly completely untouched. 

Since the filtering requirements presented by the piano 
necessitated a filter that could sweep the frequency spectrum of the 
piano, it was decided to investigate this area of tuneable digital 
filtering. Although the remainder of this naner deals with the design 
of the proposed tuner, the stress was laid on the theory and imnle- 


mentation of this filter. 


B01 aying the bearing" refers to the tuning of one octave, usually 
the one below A4 (i.e., C3-C4), by measuring its fundamentals very 
precisely and then using their upper partials for all subsequent 
tuning. This usually results in less than perfect results. See 
Refs. MH-5, MW-22 and MW-25. 
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II. DIGITAL FILTERING 


Digital filtering is the process of spectrum shaping using items 
of digital hardware as the basic building blocks. Thus the aims of 
digital fi leeting are the same as those of continuous filtering, but 
the physical realization is different. Linear continuous filter theory 
is based on the mathematics of linear differential equations; linear 
digital filter theory is based on the mathematics of linear difference 
equations. 


th 


An m~ order linear difference equation may be written as 


r 
y(nT) = ¢ L.x(nT - iT) - 


K.v(nT - iT) (3) 
1=0 7] 


nm 3 


] 


This form emphasizes the iterative nature of the difference equation; 
given the m previous values of the output y and the r + 1 most recent 


values of the input x, the new output may be computed from (3). 
Physically, the input numbers are samples of a continuous waveform 
and real-time digital filtering consists of performing the iteration 
of (3) for each arrival of a new input sample. Design of the filter 
consists of finding the constants K. and L; to fulfill a given filter- 
ing requirement. Real-time filtering implies that the execution time 
of the "computer program" for computing the right side of (3) is 
less than T, the sampling interval. 

See Fig. 6 for a pictorial representation of (3) consisting of 
unit delays of time T, adders and multipliers. 

It will be assumed in all further calculations that the set of 
input and output samples x(nT) and y(nT) are zero for all values of 


n less than zero. 


ne 


A. THE Z-TRANSFORM 


The z-transform of a sequence x(nT) is defined as 


satnt)z™ (4) 


MUL) = 


um 8 


n 
where ia is the unit delay operator defined in the s domain as 
js eee (5) 
where S = o + ju. 
For many sequences, the infinite sum of (4) can be expressed in 


closed form. For example, the z-transform of the sequence 


xinl} = On fer n<0 (6) 
x(nT) = 1. for n>0 (7) 
is 
(ee (8) 
n=0 


The transform variable Z is, in general, a comnlex variable and 
X(Z) is therefore a function of a complex variable. 


The transfer function of the filter may now be written as 


H(Z) = Z(output)/Z(input) (9) 
m Ys 

= KZ ee eee a (10) 
n=0 n=1 


The coefficients of the Z~" terms correspond to the value of the 
weighting sequence at t = nT, where n is an integer. This transfer 
function, in order to be physically realizable, must not contain 

any positive power in Z. A positive power would indicate a prediction 


or simply that the output signal precedes the input. This condition 
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implies that mr. When m =r, as in this project, a must not be 
zero. In order to ensure this condition is met, L. has been set equal 
to one in (10). 

As asserted before, the problem is now one of finding the proper 


coefficients La and K 


B. CALCULATION OF THE WEIGHTING COEFFICIENTS 


In the z-domain the zeroes may be written as 


_ JuT 
Zag = Eee (11) 
and the poles as 
JwT 
Z. Roe (12) 
and the recursive filter transfer function . then be written as: 
(WZ 92 (1-Zyg2"')---(1-22 77) 


wal 


ce Ee pel 
(1 = 12 ) (l- oy of )---01 “t, me ) 


As stated above both r and m have been set equal to n in order to 
meet the conditions for realizability. 

There remains now the step of multiplying out (13) and matching 
the derived coefficients with those of (10). The desired coeffi- 


cients are: 


Ky = 45) Pap * Sgeet --- * 26h (14) 

K, = sum of the products of the oa (15) 
taken k at a time. 

Ky ~ £51 4024037-~ “on (16) 

and 

Ce Z0 + ) + Z,3 + --- Zon (17) 

L, = sum of the products of the Zon (18) 
taken k at a time. 

Ly ~ £41 2024037-~“on (19) 


2] 


Due to the binomial coefficient nature of these terms and due to 
the fact that all of the terms are complex variables, a Fortran pro- 
gram was written to determine the various L's and K's for n as high 
as 50. Both the program and the calculated coefficients may be found 
in the computer section of this work. 

The frequency response of the designed filter may be obtained 


from: 
Hie) =) Kee! 7 (14 : Lee) (20) 
n=0 n=0 
This equation is also implemented in the program mentioned above. 
In selecting the required poles and zeroes, two apparent limi- 
tations had to be imposed: 
1. The poles and zeroes had to be real or complex conjugates 


in order that the desired coefficients would have real 
values. 


2. The poles had to lie within the unit circle to vroduce 
peoraple f1 iter. 


From the theory of sampled-data systems [Refs. EK-20, EM-27 and ER-34] 
it can be shown that the nole-zero configuration of Fig. 7, when imnle- 
mented, results in a passband with a center frequency determined 
by 
UP = Vee / (2nT) (21) 

The first factor to be determined is the distance at which to 
place the poles and zeroes around the center angle 0. in order to 
obtain the desired response. A complete investigation of this 


problem was made by Mooney in Ref. EM-28; the results will only be 


briefly summarized here. 


Ge 


Mooney found that close spacing of the poles and zeroes produced 
a sharp attenuation outside of the passband. This may easily be visua- 
lized if one thinks back on the effect the placement of poles and 
zeroes in the s-domain has on the sharpness of a continuous filter. It 
was also shown by Mooney that the passband was relatively insensitive 
to the location of the zeroes between the origin and the pole. How- 
ever, pole placement was critical and for the proper response should 
be placed as illustrated in Fig. 7. For graphs and drawings of how 
different placement affects the response of the designed filter see 


the cited reference. 


C. A TUNEABLE FILTER 
This author reasoned that since 
6 = wl (22) 
than it should immediately follow that 
Otel wT (23) 


where iy and We are the respective lower and upper cut-off frequencies 


of the passband. Subtracting (23) from (24) gave 


8 - 8 = (wo - wy )T (25) 
Ae = Awl (26) 
Ay = Sern (27) 


where Aw could now be defined as the bandwidth of the filter. 
It.was then assumed that Ae could be fixed, that that the voles 


and zeroes could be "hard-wired" in place in the z-domain. 


a3 


Equation (27) would then reduce to 
Au = C/T (28) 
or 


where fig is the sampling frequency, C and C, are constants determined 
by ae and BW is the bandwidth of the filter. 
Following the same reasoning it was easily shown that 


es wel (30) 


6o0 _ 
—— 360 Lis ij Li (31) 


where fF is the center frequency of the filter and 6. the associated 
angle in the z-domain. Once again making the angles 84 and 85 
fixed in the z-domain, (31) could be written 


ue = CA/T = C3f. (32) 


Equations (29) and (32) are the ones of prime interest. Equation 
(32) illustrates that the center frequency f. is directly related 
to the sampling frequency f.. Therefore by changing 2 it should be 
able to cause the filter to sweep the entire frequency spectrum with 
a bandwidth determined by (29). However, (29) immediately illustrates 
the existence of the ever-present "trade-offs" encountered in design. 
That is as the filter is caused to sweep and tune to ever increasing 
values of the frequency spectrum, the bandwidth is increased. Although 
the "creen" is relatively small for small values of Ae, it must be 
considered. 

This creeping bandwidth may be controlled in several ways. If 
the bandwidth requirements of the filter are not too stringent or 


the "0" of the "circuit" is not required to be too larae, this 
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widening may be neglected or taken into account when designing. This 
was the case in this project. Even though the uppver partials of the 
piano are not integer multiples of the fundamental, they are some- 
where in the vicinity of the pure harmonic frequency and a relatively 
"relaxed-0" filter could be implemented. 

Using this reasoning as a basis, Ae was selected such that as the 
sampling frequency was increased, the filter would continue to nass 
only the one desired partial of the complex sound generated. This 
is well illustrated in the computer output from program two. 

However, if the filtering requirements and "0" are such that no 
amount of variance in the bandwidth can be tolerated, the problem 
becomes more complex and much more interesting. Fiq. 8a illustrates 
two tuneable digital filters in cascade and separated by an inverse 
digital filter. The inverse filter merely changes the output of 
filter #1 into an analog form in order that the resultant wave may be 
sampled by filter #2. The shaded area of Fig. 8b depicts the spectral 
output of the second filter. It would be desireable to keep this 
interval constant as it was swept up and down the frequency spectrum. 


Let 6 be a constant and represent this interval. Then 


fies = 6 (33) 

fo = fo - 6 (34) 
but from previous discussion 

Tame inad .9 pg; (35) 

f, 85 fo) / 2n (36) 


substituting (37) and (36) into (35) 


26 


03 foo [ea Set / 2n - 6 7") 
Li, = nr / 82 (65 f.,/en-6) (38) 


which can be reduced to 


f 60 mf - b (39) 
where 

mM = 65/83 (40) 

b = 2ns/0, (41) 


These equations illustrate quite clearly that a linear relationship 
between Lis and foo can be established that results in a definite 
and invariant bandwidth. To calculate the center frequency nn of 


this passband it is clear from Fig. 8b that 


ue = fo - 6/2 (42) 
65 f.4/én - 8/2 (43) 
which results in 
un 7 af ey ad (44) 
where 
a = 0/2 (45) 
_ * 
d = “/2 (46) 


Equation (44) like (39) clearly demonstrates the fact that the center 
frequency of the passband is a linear function of the sampling fre- 
quency Foye Therefore it is quite possible to implement a tuneable 
digital filter that has an invariant bandwidth by simply adding 
another filter with a different sampling frequency which is linearly 


related to fore 
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The computer programs in the back of this paner were written to 
test the above theory of a possible tuneable fixed-coefficient digital 
filter. The first one takes as inputs the locations of the poles and 
zeroes in the complex z-domain and calculates the proper filter coeffi- 
cients. This can be done for filters up to and including n = 50. 

Filter responses for five different sampling intervals were then plotted 
via the DRAW sub-routine. The second program calculates the samnling 
frequency required for filtering the desired partials of the fundamen- 
tal frequencies of the piano. The key number, the fundamental fre- 
quency, the partial to be filtered, the maximum allowable bandwidth 

to filter the partial, theta, delta theta, the required sampling 
frequency and the resultant bandwidth of the filter are all calcu- 
lated, outputed and tabulated via this program. 

Not only does this output clearly illustrate the procedure in- 
volved in tuning the piano, but it also gives a good insight into the 
tuneable filter. That is, it demonstrates that once the coefficients 
of the filter have been set (i.e., theta and delta theta), the filter 
can be made to sweep the frequency spectrum by changing the sampling 
interval. 

The output from ovrogram number one, however, offers the most con- 
clusive proof that a tuneable digital filter is realizable. As ex- 
plained previously the program calculates the coefficients for an nth 
order filter. It then takes these coefficients and using them in the 
equations from sampled-data theory, computes the frequency response of 
the filter for five different sampling frequencies. The graphs on 
pages 46- 51 are self explanatory and offer convincing evidence that 
the previously derived equations and theory for a tuneable digital 


filter were correct. 
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D. ERRORS THAT MUST BE CONSIDERED 

Since in the process of tuning the designed filter, the sampling 
interval was allowed to move over quite a large range, care had to 
be taken in avoiding the well known “aliasing” or "foldover" distor- 
tion problem. This distortion results when the samnling frequency 
-is not at least twice as high as the highest frequency comnonent found 
in the total input signal. This problem can be overcome by employing 
wide-band samoled data filters and "pre-warping” the frequency being 
filtered [Ref. EG-15 and EW-48].° This approach was not used because 
of a much easier solution that resulted from the restrictive parameters 
of the problem. 

The highest frequency from the piano that was required to be fil- 
tered and measured was the fundamental of C8 (4186.009 Hz). Two nre- 
cautions were taken to insure that foldover distortion did not occur. 
The circuitry involved with the microphone input was designed to act 
as a low-pass filter with the upper cut-off frequency just above the 
C8 fundamental frequency. Secondly, in the design of the filter, 
be was set equal to six times the value of the highest input frequency 
possible from the "low-pass filter." 8, was then calculated from 
-(31) and "hard-wired" into the z-domain. This resulted in a safety 
factor of at least five over the entire frequency spectrum to be 


measured. 


apy wide-band sampled-data filter is meant a filter whose fre- 
quency range approaches half of the samoling frequency. The cited 
reference gives an excellent account of this type of filter and 
its possible implementation. 
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The angle 6. was computed as being equal to 30°. ae was made 
equal to 9° in order to make the bandwidth of the filter as wide as 
possible and yet still keep it within the maximum allowable range. 

Ten coefficients were called for in order that the filter would have 
as sharp a response as possible but yet still steer clear of a delay 
and coefficient accuracy problem caused by the use of too many coeffi- 
cients. The reader is referred to the computer outputs previously 
mentioned for the tabulation of all coefficients. With the mention 

of coefficient accuracy it seems avpronriate to say something about 
round-off errors, word length requirements, coefficient accuracy and 
associated problems. 

When a digital filter is realized with a digital arithmetic 
element, as is the case of the proposed tuner under study, additional 
considerations are necessary to describe the performance of the filter. 
There are three obvious degradations: 


1. quantization of the input, 


2. quantization of the coefficients of the difference 
equation, 


3. quantization of the results of the comnutations. 

All three types have been thoroughly investigated in Refs. EK-21, EM-27, 
ET-44, EK-19 and ET-45, and only the various results that affect this 
problem will be presented here. 

The two prime considerations that enter into the selection of the 
input quantization size, q, are a minimum detectable level of the sig- 
nal, Xrus and a saturation level, Xonr From Ref. EW-48 it was found 
that 


eee Sia ae 
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The A/D converter is then required to have a minimum accuracy of N 
bits, which is determined by: 
age lhe (48) 
N = log, (Xon7/Xqy) - (49) 


One word time is then defined by N data-bit times plus one sign-bit 
time. The speed of the converter, although no problem in this pro- 
posed design, may be calculated from the following formula where 

T = sampling interval 


Suz (2/Tde idee 10g5(Xen7/X7y))- (50) 


The third quantization problem mentioned is usually set equal to 
the quantization of the innut and was done so in this design. See 
Ref. EW-48 for further amplification. 

The second type of quantization error listed anpears to be the 
most important one to this writer. It is quite easy to understand 
and a thorough investigation is conducted by Knowles and Olyacato in 
Ref. EK-19. This quantization of the coefficients changes them slight- 
ly, resulting in a new, slightly different filter response. This 
happens only once when the filter is first designed. One important 
aspect of this problem is that the more coefficients that are used, 
the more accurate they must be and therefore the more this error enters 
into the problem. This is why a filter of the order of only n = 10 
was chosen. Kaiser, in Ref. EK-21, also brings out the fact that as 
the sampling frequency is raised the coefficients must become closer 
and closer to the "ideal" ones; that is, ones that would require an 


"infinite" number of bits for their implementation and storage. This 
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was the reason that the sampling frequency was set to be only six 
times the value of the highest input frequency instead of ten to twelve 
times as is done in most sampled-data networks. 

It should be pointed out here that, although the sampling 
frequency was kept as low as possible and the order of the filter 
was kept relatively small, a great deal of difficulty was encountered 
with the coefficient accuracy problem. 

The coefficients given in the computer output section were calcu- 
lated using complex double-precision arithmetic on the IBM/360/67. 
These same numbers were obtained by two completely different algorithms 
so there is little doubt that they are correct. However, when the coef- 
ficients were used to obtain the frequency response of the desired 
filter, the undesirable effects shown in Graphs 1 and 2 in Anpendix A 
were obtained. This was, however, not necessarily due to inaccurate 
coefficients. Two possible inaccuracies could have been generated by 
the nature of the library sub-programs used. First, in order to take 
the real part of the complex coefficients, the precision had to be re- 
duced from double to single. This naturally involved some tyne of 
round-off error. This necessitated change in precision resulted in 
the use of the single=precision sine and cosine library tables and 
could have introduced another round-off error. Short of actual imole- - 
mentation of these coefficients in a digital filter, any further 
attempt to obtain the frequence response by this method was imoossible. 

In order to bypass this problem Equation (13) was programmed. 

As is readily apparent by referring back to (13), this equation 
greatly reduces the probability of round-off errors. This reduction 


is caused by the knowledge of the precise locations of the noles 


3] 


and zeroes and the requirement of having to introduce only one phasor 
calculation. This is unlike the previous example where the order of 
the Filter determined how many different phasors had to be solved. 
The results illustrated in Graphs 3-6 in Appendix A conform to all 
theory that has been presented. With this discussion, this author 
feels that if the tabulated coefficients were actually implemented 


in a digital filter the frequency response of Granhs 3-6 would be 


realized. 
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IIT. A PROPOSED DESIGN FOR A DIGITAL PIANO TUNER 


A. INTRODUCTION AND OVER-ALL VIEW OF THE MACHINE 

The third and final section of this naper is devoted to the 
actual implementation of the proposed tuner. Fig. 9 illustrates 
the general flow of signals and information in block diagram form. 


ql is selected on the front 


In actual operation the key to be tune 
of the control panel by the rotation of a 12-position switch repre- 
senting the twelve notes of one of the equal-tempered octaves. 
Another rotary switch is used to select the anpropriate frequency out- 
put of the Reference Frequency Generator Module (RFGM) as will be 
explained later. The RFGM is nothing more than a very stable, crys- 
tal-controlled I.C. pulse-train generator. The output frequency of 
this module had to be at least four times as high as the desired fre- 
quency to be measured. This was required in order that the intervals 
generated by the ring counter would be quarter neriods and the subse- 
quent error count could therefore be read in cycles. 

The Time-Base Generator Module (TBGM) was added to provide the 
error count readout directly in hertz. This module acts as a real—time 
clock (RTC) and was provided with the capability of producing timing 


pulses every 5 or 10 seconds. 


Onctually the use of the word "key" here is a misnomer. AIl keys 
have more than one string in order to increase their output volume. The 
normal modern-day piano has, starting at the upper octaves and working 
towards the extreme bass, 60 keys with three strings, 18 keys with two 
strings and 10 keys with one string. Some of the larger concert pianos 
will have more but all pianos have the same 88 keys. The different strings 
are each tuned separately. This is achieved by inserting a tuning wedge 
between the strings in such a way that only one is allowed to vibrate. 
Extreme care must be taken that the multiple-stringed keys are tuned 
correctly if undesirable "beats" are to be avoided. 
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The Sign-Detection and Pulse-Shaping Module tests the left-most 
bit of the output register of the digital filter and when a positive- 
going change occurs, that is, when the bit changes from a "OQ" to a 
"1" a "one-shot" is triggered. Since an output pulse from the one-shot 
occurs only when the filtered partial to be measured crosses the 
"Zero axis,” the pulse is in synchronization with the unknown frequency 
of the filtered partial. This pulse from the one-shot is fed into 
the Frequency Deviation Detector and Counter Module (FDDCM). 

The Frequency Deviation Detector and Counter Module is a logic 
circuit that is capable of measuring the beat frequency between the 
reference frequency and the unknown input. This beat frequency or 
error count is displayed on the control panel by a series of digital 
display lights. The FDDCM also determines whether or not the beat 
frequency is siete (fast) or flat (slow) with respect to the reference. 
This enables the tuner to know in which direction to turn the tuning 
pin in the piano. For the laymen, this sharp-flat meter would be 
all’ that would be required for successful tuning. However, most pro- 
fessional tuners like to sharpen or "stretcn" the upper octaves by 
about five cents per octave to add a certain amount of "brightness" 
to the sound generated [Ref. MH-8]. This requires some form of exact 
frequency measurement, which the FDDCM does to a very precise degree. 

As mentioned previously, the input microphone and associated cir- 
cuitry was designed to act as a low-pass analog filter with an upper 
cut-off point at approximately 5 KHz. 

The Digital Filter Module multiplies the sampled input by the appro- 
priate coefficients. This involves the implementation of shift regis- 


ters and some associated logic. 
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The Display and Control Module (DCM) houses the various selection 
Switches and lights for display. Four 10-counters were properly > 
implemented and decoded to display the error in the incoming frequency 


from the piano. 


B. THE POWER SUPPLY 

The power supply had to be designed to provide more than an 
ampere of current at 3.6 volts DC to the various logic chins. This 
had to be supplied with a very low ripple factor in order to nrevent 
confusion between the ripple and the actual timing pulses emoloved in 
the circuit. It had to also provide other lower-current aiigeittiats for 
additional circuitry at +6, -6, and +12 volts DC. In addition to 
the various bias requirements it had to meet, it was also desirable to 
produce a 60-Hz signal for the Time-Base Generator Module and an AC 
voltage to the four digital-display lights. 

In order to meet all of these requirements from a single trans- 
former, eleven diodes had to be used. See Fig. 10 for the schematic 
of this module. 

The +12-volt supply was obtained from a voltage doubler consisting 
of Dl, d2, Cl and C2. The full-wave rectifier made up of D3, D4 and 
C5 provided the -6 volts, while a second full-wave rectifier, consisting 
of D5 and D6 fulfilled the +6 volts bias requirement. The +6 was 
reduced by 09, DIO and D1l1 to provide the +3.6 volts required by the 


Wee a Chins . 


C. THE ANALOG TO DIGITAL CONVERTER 
The A/D module is illustrated in block diagram form in Fig. 11. 


The operation of the ramp converter is initiated by means of a samnle 
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gate pulse which is applied to the output of the ramp generator. The 
ramp generator produces the triangular waveshape shown in Fig. 12, wave- 
form A. As this rising voltage becomes equal to or greater than the 
zero reference applied to comparator No. 1, a nulse is produced which 

is used to set a control flip-flop. When the rising voltage becomes 
,equal to or greater than the analog voltage applied to comparator No. 2, 
a pulse is produced to reset the control flip-flop. The result of this 
operation is the production of a control pulse, the duration of which 

is proportional to the amolitude of the input analog siqnal. This con- 
trol pulse is then used to gate clock pulses (which are a periodic 

train of timing pulses) into a counter which starts at a count of zero. 
The number of pulses allowed to enter the counter is proportional to 
the pulse width of the control pulse. Therefore, the final content of 
the counter is proportional to the amplitude of the analog signal. 

After the ramp voltage reaches a maximum, it is returned to its initial 
value. The trailing edge of the ramp is used to generate a reset pulse 
which reads out the contents of the counter into an output register 

and resets the counter to zero. All of the associated waveforms are 
illustrated in Fig. 12. The "transfer" occurs when the seven-bit 
counter is read-out into the read-out gates. 

The seven-bit counter illustrated in Fig. 11 had in reality eight- 
bits. The left-most bit, the sign-bit, was dropped from the above 
discussion for the sake of clarity. The seven-bit counter (including 
the eighth-sign-bit) was implemented because it produced a very high 
and desirable degree of accuracy. That is, it was capable of dividing 


7 


the input signal into 628 levels (2° = 628). 
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D. THE FREQUENCY DEVIATION DETECTOR AND COUNTER MODULE 

If two frequencies are superimposed on an oscilloscone screen, 
the relation between them can be determined by counting the number of 
times they move in and out of synchronism. In reality this nothing 
more than determines a beat frequency, which will give the deviation of 
the second from the first, if the first is considered to be the refer- 
ence. In existing analog systems the attempt to realize this beat 
frequency at low frequencies and with small errors can be quite complex 
and cumbersome. For example, an error of 0.01% at 400 Hz would produce 
a beat of one twenty-fifth of a cycle per second. Although accuracy 
this high is not demanded in the tuning of the piano, a reasonable degree 
of accuracy is required and can be realized with a saving in size, 
cost and required skills of an operator by the use of digital circuits. 
Refer to Fig. 13 for the logic diagram of the designed FDDCM. 

The crystal-controlled reference frequency from the RFGM is fed 
into the phase generator, a simple four—stage ring counter (See Fig. 14), 
to produce the required timing intervals. Four intervals, A thru D, 
were provided to allow for the polarity of the errors to be calculated. 
The one-shot (mono-stable multivibrator) produces a nulse that is in 
Synchronization, as explained previously, with the unknown frequency 
to be measured. This pulse, designated fy was ANDed separately with 
each of the four phases from the ring counter. This gives the outout 
of the coincidence detector as (AF, BF, CF and DF). These signals 
provided the basis for the remaining logic. The error-count detector 
was an RS flip-flop which produced an error count whenever the unknown 


had cycled from coincidence with phase A to nhase C. The recordina of 
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an error count occurred only when the output of this RS flip-flop 
changed from "0" to the "1" state. Note that repeated pulses on the 
set input can produce no further error counts. 

The polarity, that is whether or not the unknown is fast or slow 
compared to the reference, was determined by the remaining logic. 
Fiqs. 15, 16 and 17 illustrate the various waveshapes for a "sharp," 
"flat" and perfectly tuned unknown input. If the input is perfectly 
"tuned" to the reference it will appear to remain stationary when 
plotted against the reference with respect to time. However, if 
jt were sharp or flat, it would tend to drift to the right or left, 
respectively, and thus cause an error count to be generated. 

Once again, only off-the-shelf Fairchild Semi-Conductor Micro- 


Logic Chips were employed in all design work. 


E. THE REFERENCE FREQUENCY GENERATOR MODULE 

The RFGM works on the same basic principle as the pitch reference 
found in Ref. EL-22. That is, it takes a set, crystal-controlled 
frequency and uses it to drive a chain of JK flip-flops. This chain 
divides the crystal frequency output into the desired reference fre- 
quency outputs. 

The power supply used in the cited reference was deleted in favor 
of the one already discussed. This was only for the sake of simplicity 
and in an attempt to avoid needless repetition of components. The 
Speaker and associated circuitry were also dropped for the same reasons. 

The basic block diagram of the circuit appears in Fig. 18. When 
one of the crystal oscillators is switched on, the output wave is taken 


by the input of the Logical Schmitt Trigger, illustrated in Fig. 23, 
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and made into a square wave of the same frequency. Take for example, 
the switch Sy as being in the position indicated in Fig. 19. Then 

the output of the Schmitt Trigger is a 34.29 MHz square wave and is fed 
into the chain of JK flip-flops. With the switch in this position there 
are 12 flip-flops in this chain when tuning octaves five, six and seven. 


'2°_ 4096). The 


With no feedback, the chain will divide by 4096 (2 
resulting output frequency of the chain would therefore be 8371.9 Hz. 
Notice that this frequency is four times the frequency of the desired 
fundamental of C7, the desired first partial of C6 and the third partial 
of C5. It is therefore the proper reference frequency required to 
input to the FDDCM when tuning these notes. 

Now if the frequency output from the previous example were buffered 
and fed back into the outputs of flip-flops one, two, five and six, 
the chain would now divide by 3866. This is because extra counts 
exactly equal to the difference between the desired divisor of 3866 
and 4096 have been added. In arithmetic form the preceeding is 
illustrated by 

4096 - 3866 = 130 
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When this division by 3866 is carried out on the input frequency of 34.29 
MHz, the resulting output was calculated to be 8869.6 Hz. This is four 
times the desired fundamental of C7#, the desired first partial of C6# 


and the desired third partial of C5#. Once again the required 


oie, 


reference frequency to tune these notes has been obtained. This same 


procedure was carried out for the 36 notes in the upper three octaves 


and the results are displayed in Table 2, | 
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NOTE DIVISION DESIRED PARTIAL REFERENCE 
RATIO FREQUENCY FREQUENCY 
C7-C6-C5 4096 2093.0 8371.9 
C7#-C6#-C5# 3866 DVB 8870.0 
D7-D6-D5 3650 2349 .3 9397.2 
07#-D64-D5# 3444 2489.0 9956.0 
E7-E6-E5 3250 2637.0 10548 .0 
F7-F6-F5 3068 2793.8 11175.2 
F7#-F64-F5# 2896 2959.9 11839.6 
G7-G6-G5 2734 3135.9 12543.6 
G7#-G64-G54 2580 3322.4 13289.6 
A7-A6-A5 2436 3520.0 14040.0 
A7#-A6#-A5# 2298 3729.3 14917.2 
B7-B6-B5 2170 3951.1 15804.4 
TABLE 2 


Reference frequency output data for 
the upper three octaves 


As is illustrated in this table all of the upper three octaves of the 
piano can be tuned with the 34.29 MHz frequency setting and the chain 
of twelve dividing flip-flops. 
As explained previously it was desired to correctly measure and 
tune the fourth partials of the notes in the fourth octave. This octave 


of frequencies can be obtained by dividing the third nartials of the 


"Vai though the reference frequency output is not exactly four 
times the value of the partial being measured, it is accurate to 
+ 0.5 cent, making it twice as good as the best tuning fork 
available. 
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fifth octave by two. Since it has already been shown how the RFGM is 
able to generate the correct frequency for the fifth octave, the only 
requirement is to divide this set of frequencies by two in order to 
obtain the correct set of frequencies for the fourth octave. That is, 
when tuning the fourth octave, another flip-flop must be switched in 
from the control panel. This is done by simply rotating the seven- 
position switch marked "octave" to the 4th position. 

It was also desired to measure and tune the fourth partials of 
the lower three and one-third octaves. It would have been highly 
desirable to simply continue adding divide-by-two flip-flops in the 
dividing chain to decrease the reference frequency to its desired 
value. Had this been done, however, it would have been required to 
filter and measure the ninth nartials of the third octave, the eighth 
partials of the second octave and the seventh partials of the first 
octave. According to Plomb's criterion this would have presented no 
problem. However, these specific partials are unfortunately sup- 
pressed by manufacturers because of the shrillness they tend to add 
to the music when present [Ref. EW-221.'¢ Clearly, another approach 
was indicated. 

This observation required that a second crystal with a resonant 
frequency of 6.7171 MHz be added. When this frequency was divided by 
the basic chain of twelve flip-flops, a set of frequencies was 
produced that was four times the desired partial frequencies of the 


fourth octave. Therefore, the required reference frequency had once 


Ler is done by having the hammer hit the string between one- 
seventh and one-ninth of the distance from the one end of the string 
to the other. 
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again been obtained. The first and second octaves can also be tuned 
from this crystal by switching one more flip-flop for the second oc- 
tave and two more for the first. Once again this is achieved by the 


rotary switch on the face of the Display and Control Module. 


F. THE TIME BASE GENERATOR MODULE 

The TBGM, illustrated in Fig. 19, consists of seven JK flip-flops 
and one dual two-input gate package (Fairchild Semiconductor RTuL 9914). 
Incoming 60-Hz clock pulses from the power supply are fed to FFl and 
FF2 which divides the input by three. The resulting 20-Hz signal 
provides clock pulses at 0.05-second intervals (T = 1/f). FF3 divides 
the 20-Hz signal by two to provide clock pulses at 0.10-sec intervals, 
and FF4, FF5 and FF6 are connected in a divide-by-five circuit config- 
uration to obtain the 0.50-second timing interval. The 2-Hz output 
of FF6 is divided by FF7 for the 1.0-second clock pnulse. The first 
half of the dual two-input gate is used as an inverting amplifier to 
boost the divide-by-five input signal level and the other half is 
used to amplify the final clock—pulse output. This outnut chain of 
pulses is fed into a divide-by-ten counter illustrated in Fig. 20. 

The divide-by-ten module was designed using the set and clear of 
the basic JK flip-flop as gates. It was able to feedback directly to 
the binary divider without any extra parts and inhibit counts 11 
through 16. Although the circuit is quite difficult to decode, it 
provides a very simple, stable and cheap divide-by-ten package. 

The output from this circuit provides the actual five or ten-second 
timing pulse. When the pulse occurs it gates the FDDCM and halts all 


counting. 
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In order to ensure that the 60-Hz signal obtained from the power 
supply had a relatively sharp rise and fall time, a Logical Schmitt 
Trigger was implemented between the power supply and the input to the 
TBGM. It is illustrated in Fig. 21 and as shown emnloys a Fairchild 
RTyL 9914 gate. The series resistance limits the peak voltage at the 
gate input to two to three and one-half volts. The capacitor acts 
as a filter and the diode clamps the negative portion of the wave to 


ground. 


G. THE TUNEABLE DIGITAL FILTER MODULE 

The implementation of this module might have presented a real 
problem if the use of a tuneable filter had not proven feasible. It 
would have then been necessary to generate a new set of coefficients 
for each octave to be tuned. Fortunately, as has already been dis- 
cussed in some detail, this was avoided. 

The filter used has the same characteristics as the one discussed 
in Ref. ET-43, except that the sampling time could be varied by the 
tuner from the CDM. This Fixed-Coefficient Diode-Array Digital Filter 
was employed because of its extremely small size and simplicity of 
implementation. It uses the table-lookup method of product generation 
where the diode-array is used to store the multiplication table. A 
diode-array mechanization is illustrated in Fig. 23. 

The input data from the A/D converter is shifted in serially to 
load the shift register. The contents of this register addresses the 
product, which is stored, in the diode array. The product is loaded 
into the output reqister which supplies the input to the next chip or 


array. Each chip delays the information the required interval T. 
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The schematic of Fig. 22 can be implemented on a single MOS chip. 
For this proposed tuner it would take 20 of these chins; one for each 
numerator and denominator coefficient. Reference ES-38 discusses 
recent developments in this field and also the possibility of putting 


many of these arrays on a single chip. 


H. THE DISPLAY AND CONTROL MODULE 

The DCM consists of a set of digital display lights, a meter indi- 
cating flat, zero and sharp, and two sets of selection switches. The 
selection switches control the number of flip-flops in the RFGM dividing 
chain and the sampling interval of the A/D. The meter is a three— 
terminal galvonometer and will tell the tuner whether the partial 
being measured is sharp or flat with respect to the reference fre- 
quency. The display lights will measure a frequency deviation up to 
999.9 Hz. There are four lights on the panel and each one is used 
to decode a "10" counter. The counter and the required decoding 
scheme are illustrated in Fig. 24. The output from the first counter 
is used as the trigger for the second. The output from the second is 
used as the input for the third and the output from the third is used 
as the input for the fourth and final stage. Thus, a "1000" counter 
has been effectively achieved. When employing the usual 10 second— 
timing interval, the counter becomes a "100" counter that counts by 


tenths. 


fe 


IV. CONCLUSION 


From the study of the physics of the piano, it has been shown 
why it was necessary to measure and tune the eigen-partials of the 
complex sound produced by the steel-stringed piano. It was also 
demonstrated how this filtering and frequency measuring could be 
achieved with a tuneable fixed-coefficient digital filter and several 
additional I.C. modules. The remaining portion of the work was then 
devoted to the possible implementation and construction of these 
modules. Based on the sound and music theory presented this pro- 
posed design seemed to offer the best electronic tuner to date. 
Unfortunately time did not allow for the actual interfacing of the 
modules and testing of the tuner as a single unit. 

However, it should be stated that this author in no way advocates 
this machine as a panacea for the complex problem of tuning equal- 
tempered key-board instruments. Music and its desirable qualities 
are things which surpass scientific measurement because, contrary to 
all attempts to define them otherwise, they remain purely subjective 
quantities. A piano tuner who employs his well-trained ear is able 
to take the tastes of the period into account and tune accordingly, 
This machine or any other can only give him a better or more per- 


ceptive "ear." 
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Graph 2. Illustrating the Results of Coefficient 
Inaccuracy on the Frequency Response of a 
Digital Filter. 
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Graph 3. Frequency Response of the Designed 
Digital Filter, 
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Graph 6. Frequency Response of the Designed 
Digital Filter, 
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Fiaq. 2. Vibrations of an Ideal String. 
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